Mixing of ultracold atomic clouds by merging of two magnetic traps 
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We demonstrate a method to make mixtures of ultracold atoms that does not make use of a two- 
species magneto-optical trap. We prepare two clouds of ^'^Rb atoms in distinct magnetic quadrupole 
traps and mix the two clouds by merging the traps. For correctly chosen parameters the mixing 
can be done essentially without loss of atoms and with only minor heating. The basic features of 
the process can be accounted for by a classical simulation of particle trajectories. Such calculations 
indicate that mixing of different mass species is also feasible, opening the way for using the method 
as a starting point for making quantum gas mixtures. 
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I. INTRODUCTION 

Production of ultracold atomic gases consisting of a 
mixture of two different elements is interesting for several 
reasons. First of all it makes it possible experimentally to 
study and exploit inter-species interactions [l|, [3, H, 0, H, 
HSJa HJlO] including inter-species Fcshbach resonances 
|lllll2l.ll3i,,14]. An important possibility in this respect is 
to produce ultracold heteronuclear molecules with strong 
electric dipole-dipole interactions using photoassociation 

fi, 16, 17, 18, 19] or Fcshbach resonances ^W, 20, 21, 2l, 
. Another is the study of phase separation phenomena 
HSi [Hi [IB] and realization of fermionic superfluidity 
mediated by bosons [13, [H, [11] . Furthermore mixtures 
in optical lattices is a new and relatively unexplored field 
of study with many interesting phenomena [20l [30l . [3l[ . 
[3i,[3l,[3i,[3l. 

The mixtures also have an important application in 
that one of the species can be used to sympathetically 
cool the other one. If one of the species is evaporatively 
cooled then the other species - which might have scatter- 
ing properties which make evaporative cooling impossi- 
ble - will be cooled too if the inter-species thermalization 
rate is large enough. This principle has been used to cool 
fermionic atoms to quantum degeneracy [1, 2, 3, 4, 6, 7, 8] 
and can also be used to condense bosonic isotopes which 
are otherwise difficult to condense as demonstrated for 
^iRinRef. 

Currently many groups have produced ultracold mix- 
tures in the quantum regime. This includes boson- 
fermion mixtures where the fermions are cooled sympa- 
thetically by the bosons (^Li/^Li [3, [1], ^SNa/^Li H, 
""Rb/''°K^, [1 [|, 87Rb/6Li 0) and very recentl 
He/'*He |37 ) and boson-boson mixtures (*^Rb/'*^K [3/ 
Rb/ 
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I. [l3]). Furthermore several other mixtures 
in the ultracold, but not yet quantum, regime have been 
studied (e.g. [i,[l3,[iii[l3i[ll[li,[3i). 

Even if the goal is to enter the quantum regime by 
means of evaporative cooling the starting point is a mix- 
ture which has already been cooled to the 100 /.tK regime 
by means of laser cooling. The standard approach used to 



produce such cold mixtures is to cool and trap the two 
(or even three [39]) species simultaneously in the same 
magneto-optical trap (MOT). 



In this article we demonstrate a completely new 
method for making the mixtures. The basic idea is to cool 
and trap the two elements in separate magneto-optical 
traps and mix the two clouds afterwards. After having 
cooled the atoms we trap the two clouds in distinct mag- 
netic quadrupole traps and mix them by merging the 
two traps mechanically. Since the magnetic fields from 
the two traps interfere during the merging process it is 
not trivial that one can mix atomic clouds in this way. 
We have demonstrated that the method is indeed feasible 
and a promising new way of producing mixtures of cold 
atoms, but it needs accurate control of various parame- 
ters. 



The procedure makes it possible to optimize the cool- 
ing and trapping process of each species individually and 
it is a modular approach where one can build the setup 
needed for producing one species and then add the other 
one when needed. Furthermore there is no need for optics 
which can handle two different wavelengths. 

In section |IT] we give an overview over our experimen- 
tal setup. Section IIIII is devoted to a discussion of the 
trapping potential and its evolution during the mixing 
process. In section HVl we introduce a simple theoretical 
model which we use to simulate the atomic particle tra- 
jectories in order to make predictions about the outcome 
and understand the physics of the mixing process. In sec- 
tion |V] we discuss our experimental results and compare 
them with the results of the theoretical model which turn 
out to reproduce the essential features quite well. We 
examine the dependence on essential parameters such as 
relative and absolute trap depths. Finally in section IVTl 
we discuss mixing of different elements based on the the- 
oretical simulations and conclude in section fVIII 
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Figure 1: Sketch of our vacuum chamber. The coils constitut- 
ing the magnetic traps are symbolized by the dashed circles. 
The distance from the MOT to the corner and from the corner 
to the science chamber are 49 cm and 37 cm, respectively. 



II. EXPERIMENTAL SETUP 



in the coils up to 150 A during 1 ms. Afterwards we wait 
for 55 ms and increase the current to 250 A in 200 ms be- 
fore moving the atoms through the differential pumping 
hole. 

The procedure which we use to mix the two atomic 
clouds is explained in detail in Fig. [21 When we have 
caught atoms in the traps, we move Trap 1 towards Trap 
2 at a speed of 5 cm/s. Trap 2 is stationary during the 
process. Both clouds are situated in the "low pressure" 
part of the chamber during the process. After the mixing 
process we move the atoms to the science chamber and 
transfer the atoms to the science chamber quadrupole 
trap. Here we wait 10 s to allow for thermalization be- 
fore we take absorption images after 7, 9, 11 and 13 ms 
time-of-flight to measure the temperature and number of 
atoms. 

In order to study how each of the atomic clouds is 
influenced by the process we can choose to have atoms in 
one or the other of the two traps or in both of them. We 
do this by having the cooling light either on or off during 
each of the two MOT phases. 



Our experimental setup is sketched in Fig. [TJ The 
vacuum chamber consists of two parts separated by a 
differential pumping hole: a "high pressure" part (P ~ 
2 X 10~^° torr) with a cylindrical glass cell where we cool 
and trap ^^Rb atoms in a MOT 40] and a "low pressure 
part" (P < 10~^^ torr) where the lifetime due to rest gas 
collisions is about 2 min. 

We have two magnetic quadrupole traps which are 
mounted on computer controlled mechanical positioning 
systems like in Ref. [43|. One trap is for moving the 
atoms from the MOT to the corner of the chamber and 
the other can move them from the corner to the science 
chamber (see the sketch in Fig. [T]). The science cham- 
ber is made of stainless steel and has its own magnetic 
quadrupole trap which can be turned off in 340 /is allow- 
ing us to measure temperatures by time-of-flight imaging. 
The science chamber is furthermore used to make Bose- 
Einstein condensates (BECs) by standard forced radio 
frequency evaporation in a Quadrupole-Ioffc configura- 
tion (QUIC) trap [4l|. 

The setup is prepared for the addition of a second 
source chamber at the right side of the corner and this 
part will have its own movable magnetic trap. The idea 
is to mix the two elements at the corner in the future. 

To study the principle of mixing by merging two mag- 
netic traps we mix two clouds of ^^Rb in the |P = 
1, tof = — 1 > hyperfine state. To produce the cooling 
light for the MOT we use a standard diode laser setup 
Ea, 113, El- During the last 5 ms of the MOT phase we 
compress the atoms in the MOT by detuning the cooling 
lasers and decreasing the power in the repump to mini- 
mize heating when the magnetic trap is turned on [40| . 
The trap coils provide the magnetic field of the MOT. Af- 
ter the compressed MOT phase about 40 % of the atoms 
are caught in the magnetic trap by ramping the current 



III. MAGNETIC FIELD AND TRAPPING 
POTENTIAL 

Rb atoms can be trapped using the linear Zeeman shift 
[4^ . For a magnetic field of magnitude B{r) the trapping 
potential is U{r) = g-pmFi^BB{r) where gp = —1/2 for 
the P = 1 hyperfine state and gp ~ 1/2 for the F=2 
hyperfine state of the ground term. 

Our magnetic traps are simple quadrupole traps ideally 
consisting of two circular current loops with radius R 
and separation distance 2A. Around the zero-point in the 
middle of the trap the magnetic field increases linearly 
in all directions with the axial gradient being twice the 
radial gradient. Situating the origin in the middle we will 
refer to the axis through the coil centers as the z axis, 
and the radial coordinate will be called p. Ref. [45] gives 
analytical expressions containing elliptic integrals for the 
magnetic field from circular current loops. 

Near the minimum (for our coils up to about 1 cm from 
the minimum) the magnitude of the magnetic field can 
be well approximated by [46^ 

B{p,z)^^^^Ti^ (1) 

where d\B\/dp is the radial magnetic field gradient. This 
formula shows that the potential energy isosurfaces near 
the minimum are ellipsoids and it makes it possible to 
make analytical calculations of thermodynamic quanti- 
ties f section HV Ap . 

The coils of Trap 1 and Trap 2 consist of, respectively, 
16 and 31 windings of 4.25x4.25 mm^ square, copper 
tubing with cooling water fiowing in the middle. The 
coils are cast in epoxy for stability and each coil has outer 
dimensions of about 10x10x2 cm. 
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Figure 2: The procedure used to mix two Rb clouds, (a); 
First we collect atoms in the MOT, trap them in Trap 1 and 
move the trap to the corner where the atoms are transferred 
into the other movable trap (Trap 2). This trap is moved 
7 cm towards the science chamber to shield the atoms from 
scattered light, (b): Trap 1 is moved back and a second 
cloud of atoms is collected. This cloud is moved into the 
" low pressure" part of the chamber to a position 22 cm from 
the corner. The trap current is slowly changed to the initial 
value wanted during the mixing. Trap 2 is moved back to 
the corner, (c): Finally we merge the two traps by slowly (5 
cm/s) moving Trap 1 towards the corner until the trap axes 
are 2 cm apart. This is the mixing process. We cannot move 
the traps closer together due to the limited space around the 
vacuum chamber. For the same reason Trap 2 is displaced 5 
mm from the transport axis of Trap 1. 



We have calculated the magnetic field from our coils 
using a numerical implementation of the Biot-Savart law. 
As it turns out these calculations show that for our pur- 
poses the magnetic field from the coils can be very accu- 
rately approximated by the field from two circular wind- 
ings, if A, R and / are chosen properly. The validity 
of this approximation makes the numerical simulations 
much faster and shows that if the mixing process works 
with our coils, it should work for other quadrupole traps, 
too. 

Due to the interference between the fields from the 
two traps during the mixing process, the magnetic field 
from each trap needs to be accurately determined in the 
radial direction even at distances up to 20 cm. Therefore, 
to find optimum values of A, R and J, we fit the Biot- 
Savart calculation for Bp{p) in the z = Q plane with 1 A 
through the actual coils with the corresponding function 
-Ba,r,i (p) for a single pair of circular windings using A, R 
and / as fitting parameters. This representation is very 
good: Within 2 cm from the z=0 plane the maximum 
deviation is 7 mG and this occurs in the middle of the 
trap where the numerical Biot-Savart implementation is 
inaccurate due to the discontinuity of the field gradient. 
The best fitting values of A, R and / are listed in Table 

in 



Table I: Circular coil pair parameters. 



Coil 


Trap 1 


Trap 2 


A (cm) 


3.941 


7.082 


R (cm) 


3.2 


4.21 


1(A) 


16.0 


32.17 


Radial field gradient (G/cm)/A 


0.36 


0.20 


Field barrier height (G/A) 


0.77 


0.59 



The maximum current we can use in our coils is 400 A 
which is limited by the current supply. 

We use the current ramps sketched in Fig. [31 The cur- 
rent in the coils constituting Trap 2 is constant through- 
out the mixing process. The current in the moving Trap 1 
is constant in the beginning but is linearly ramped down 
to towards the end such that it is exactly when the 
movement stops. Also shown is the position of Trap 1 
as a function of time. We use soft acceleration and de- 
celeration phases (with constant time derivative of the 
acceleration) , but during most of the time the velocity is 
constant. 

We have tried to vary the time At during which we 
ramp down the current of Trap 1 keeping the total mixing 
time constant (see Fig. [3]). We have found that the exact 
value of At is not critical as long as it is larger than 
1.5 s. Ramping down the current of Trap 1 faster than 
that causes additional heating (but essentially no loss of 
atoms). The optimal ratio of initial currents in the two 
traps is strongly dependent on the value of At. 

We have also tried to make deviations from the cur- 
rent ramp in Fig. |3|by ramping down in two linear steps 
instead of one. This, however does not lead to any im- 
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provement of the results. One could imagine that a suit- 
ably tailored soft curve would lead to marginally better 
results, but since the results using the simple linear ramp 
are quite satisfactory we have not pursued this issue fur- 
ther. 




Time (s) 



Figure 3: (Color online) Current ramps and the position of 
Trap 1 during the mixing process. 
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Figure 4; (Color online) Position of Trap 1 during the mixing 
process (blue crosses) together with the calculated position of 
the trap using the velocity and the stop time as fitting pa- 
rameters (red curve). The trap position was determined from 
digital images of the coils taken at well defined times during 
the process. The calculated position only fits the measured 
one if the stop time is correctly chosen to within about 10 ms. 



The current ramp must be accurately synchronized 
with the movement of Trap 1. We apply a reset proce- 
dure on the mechanical positioning system to eliminate 
day-to-day drifts in the trap position. After reset the un- 
certainty in the coil position is about 50 /xm. By taking 
images of the trap coils during the movement we have de- 
termined the exact position as a function of time (Fig. [4]) 
and used this information to time the current ramp such 
that the current in Trap 1 becomes zero exactly when 
the trap has reached the final position. By using the 
same procedure we can conclude that the timing drifts 
of the trap movement is less than 10 ms and this level of 
stability is indeed necessary to get reproducible results. 

Armed with a precise magnetic field calculation we are 
able to study the geometry of the trap potential during 
the mixing process. This is best viewed by plotting po- 
tential energy isosurfaces of the magnetic field. Some 
representative isosurfaces are shown in Fig. [5l The 
shown isosurfaces correspond to an energy of gprnpliX 
2 mK where k is Boltzmann's constant. Gravity is also 
included in the calculations. If the atomic clouds were 
in thermal equilibrium during the mixing process (and 
we stress that they are not) 35 % of the atoms would re- 
side outside these isosurfaces at a temperature of 200 fiK. 
Therefore it is certainly not clear from these isosurfaces 
that the mixing can be done successfully, but nevertheless 
they give an impression of what is going on. 

When the two traps are far apart there are 3 isosur- 
faces - one around the two trap centers and one between 
the two traps. The isosurface in the middle contains no 
atoms in the beginning. 

When the traps are moving closer together the poten- 
tial gets weaker and the isosurfaces begin to mix which 
makes it possible for the atoms to leak into the region 



in the middle. This is the critical point since the atoms 
are thermodynamically allowed to spread out over a large 
volume. That can eventually cause heating when the trap 
potential tightens again at the end of the mixing process 
and can lead to loss if atoms hit the wall of the 40 mm 
vacuum chamber tube. Shortly after this, all 3 isosur- 
faces are merged and in the end we are left with only one 
isosurface which ideally contains all the atoms. 

If the current Ii in Trap 1 is too small compared to 
the current I2 in Trap 2, then the isosurface of Trap 1 
mixes with the isosurface in the middle too early and the 
atoms in this trap are lost. On the other hand, if the 
ratio /1//2 is too large, the atoms in Trap 2 are lost (see 
Fig. El). This is exactly what we see experimentally, and 
fortunately there is a well defined current ratio which 
leads to successful mixing. 



IV. CLASSICAL SIMULATION OF THE 
DYNAMICS 

In order to get a theoretical understanding of the mix- 
ing process we have carried out simulations of the particle 
dynamics. The temperature of the clouds is so high that 
quantum effects can be safely neglected. One way to 
see this is to consider the typical spacing between quan- 
tum mechanical energy levels in the linear trap potential. 
From dimensional arguments, or by looking at the energy 
levels in the 1-dimensional potential V{x) = C\x\ ([47[ p. 
109) this spacing is of the order of 

/(-i2fc2\ 1/3 

E=i-^—\ =kx0.2^iK (2) 
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Figure 5: (Color online) Selection of potential energy isosurfaces during the mixing process with optimal parameters. The 
surfaces correspond to a potential energy of gFmpk x 2 mK. For F = 1 atoms at a temperature of 200 fiK 65 % of the atoms 
would stay inside these isosurfaces in thermal equilibrium (however, we stress that the atoms are not in thermal equilibrium 
during the mixing process). Notice that the time interval between the third and the fifth plot is only 0.2 s. The vertical 
asymmetry at large distances from the z—Q plane is due to gravity. The current ramps are as in Fig. [3] The current in Trap 2 
(situated at x = cm) is /2=265 A and the initial current in Trap 1 is Ji =322 A corresponding to the ratio /1//2 =1.215. 



where k is Boltzmann's constant using the worst case 
parameter C = fiB x 150 G/cm. This is 3 orders of 
magnitude smaller than the typical thermal energy. Fur- 
thermore we neglect interactions and collisions between 
the particles. See section IIV Al for a discussion of this 
approximation. 

So we propagate the particles classically: 



magnitude given by the cyclotron period 



dri 

dt 



dVi 

' dt 



-V?7(r,) 



(3) 



Ti being the position of each particle. We implement 
V/7 = g-pm-p^B^ B -I- mgz numerically and solve the 
equations using a standard multi-order Runge-Kutta 
method with error control in Matlab. The dynamics is 
rather complicated and cannot be described analytically. 
However, all particle trajectories exhibit some oscillatory 
behaviour with a characteristic period with an order of 



27r 



mp 



g-pm-p^B X d\B\ldp 



(4) 



which is about 40 ms initially. 

From the result of such a simulation we can extract 
various interesting quantities like position, extension and 
temperature (Fig. [7]) of the atomic clouds as a func- 
tion of time. It is easy to discriminate the lost particles 
from the particles which stay trapped since the former 
perform a nearly free fall in the field of gravity and will 
quickly travel far away from the trapping region. Since 
our vacuum chamber tube has a diameter of 4 cm, par- 
ticles which during the process reach a position which is 
more than 2 cm away from the transport axis are consid- 
ered lost and removed from the calculation before calcu- 
lating the quantities above. 

Using optimal parameters it turns out that if we con- 
sider the atoms originating from Trap 1 and Trap 2 in- 
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Figure 6: (Color online) Isosurfaces corresponding to a po- 
tential energy of gFrnpkx 2 mK for current ratios which are 
far from optimal. The current in the coils of Trap 2 (situated 
at a:: = cm) is I2 =265 A. The initial current in Trap 1 
(situated at a:: = 8.93 cm) is respectively Ii —280 A (upper 
figure) and 350 A (lower figure). The figure shows that if the 
current ratio 7i /I2 is too small the atoms in Trap 1 are likely 
to become lost or heated, if it is too large the atoms in Trap 2 
are likely to suffer from loss or heating. Notice that the scales 
of the axes are different from those in Fig. [S] 



dividually, then the kinetic and potential energies cor- 
respond to the same effective temperature and there is 
almost no anisotropy of the velocity distribution. This 
means that these ensembles are not far from thermal 
equilibrium at the end of the mixing process. But the 
ensemble from Trap 1 is not in thermal equilibrium with 
the ensemble from Trap 2. However, since thermalization 
takes place through elastic collisions the equilibrium tem- 
perature of the mixed cloud as a whole can be found from 
the conserved mean energy, since in thermal equilibrium 
< E >=< E^kin > + < U (9/2)fcT as described in 
section IIV Al Also the equilibrium phase space density 
can be calculated. 



A. Thermodynamics of the quadrupole trap 

Using the approximation ([T]) thermodynamic equilib- 
rium quantities of the trapped atoms can be calculated. 
We introduce the radial trap potential gradient 
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Figure 7: (Color online) Results of a simulation of the particle 
dynamics using the current ratio I1/I2 =1.215 which gives 
the optimal mixing results, (a) shows the position of the 
two clouds (the mean position of the atoms) as a function 
of time, (b) shows the RMS extension of the clouds in the 
mixing direction, (c) shows the mean energy of the two clouds 
expressed as the corresponding equilibrium temperature T = 
{2/9k) < E >. The temperature temporarily decreases just 
before the clouds merge, because the trap potential opens up. 
Since the potential energy gradient of Trap 2 is smaller than 
the initial gradient of Trap 1 the phase space density of the 
cloud in Trap 1 decreases even if the final temperature is not 
higher than the initial one. Particles which are lost during 
the mixing process have been removed from the calculation. 
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The density is given by the Boltzmann distribution: 

n(r) = £(^)%-V^/- (6) 

where N is the number of atoms. 
The peak phase space density is 



n(r = 0)Xl^ 



N f C 



4tt \kT J \2TTmkT J 



3/2 



y9/2 



The mean kinetic and potential energies are 



< Ekin >= ^kT <U >= 3kT 



(7) 



(8) 



arising solely from the fact that the energy is quadratic 
in momentum and linear in position. 

In order to determine the mean elastic collision rate 
we need to know the mean density which is given by 



< n(r) ^ y '^^^ 



3„ _ 1 



N 



n{r = 0) (9) 



The mean time between elastic collisions becomes 



Tel = (< n > (T < Wiol > 

10* / T 

^ 19 ms ■ 



(10) 

/ 1 100 G/cm^ ^ 



N \^OO^lKJ Vspmp d\B\/dp 



(11) 



using the scattering length a = 106 oq from |4^, the 
cross section a — Sttc^ and the mean relative velocity 
< I'rei >= \/ 8kT / {TTm / 2) . For the experimental param- 
eters used in this work r lies in the range 1-8 s in the 
science chamber trap. This means that thermalization 
is expected to occur after the mixing process has taken 
place. Experimentally our coldest clouds (those we get 
when we use optimal parameters) are thermalized shortly 
after the transfer to the science chamber trap while the 
hotter clouds take on the order of 5-30 s to thermalize. 

During the mixing process, where the magnetic fields 
of the two traps interfere, the mean density is smaller so 
the mean time between collisions is even larger. At a mix- 
ing speed of 5 cm/s the mixing effectively takes place in a 
fraction of a second, which implies that thermalization is 
highly unlikely to be important. Therefore it is reason- 
able to neglect collisions in the calculations during the 
mixing process. After the mixing process, however, the 
clouds begin to thermalize and experimentally we store 
the atoms for 10 s in the science chamber quadrupole trap 
before we take the image in order to be close to thermal 
equilibrium. 



B. Initial conditions 

In each trap we prepare an ensemble of 500 particles 
obeying Boltzmann statistics. We assume that the phase 



space density is conserved during transport and adiabatic 
compression so we choose the initial temperature 
such that the phase space density matches the one de- 
termined from temperature measurements in the science 
chamber. At a radial trap gradient d\B\/dp of 84 G/cm 
we measure the temperature of a cloud which has not 
gone through the mixing process to be 195 fiK. Using 
the scaling of the phase space density ([7]) we see that the 
proper initial temperature is 



84 G/cm 
d\B\/dp 



-2/3 



X 195 fiK 



(12) 



so the two clouds have different temperatures but the 
same phase space density. 

The initial velocities are found by selecting each Carte- 
sian component randomly from the Gaussian distribu- 
tion: 



P{vj) cx e' 



J = x,y,z 



(13) 



We distribute the positions using a rejection method, see 
Footnote 1501. 



V. RESULTS 

Using optimal parameters we get the following main 
result: If we mix two clouds with about 1.7 x 10* atoms 
each, we loose only about 5% of the atoms and the mea- 
sured temperature in the science chamber increases from 
195 fjK to 310 fjK. This means that according to Eq. 
([7|) the phase space density relative to the phase space 
density of the initial clouds is 24 %. Using optimal pa- 
rameters in the theoretical model we get a relative phase 
space density of 27 % in excellent agreement with the 
experimental result. The elastic collision rate is also im- 
portant in relation to subsequent evaporative cooling. It 
scales as TV/T^/^ (see Eq. ^) which means that the 
mixed cloud has an elastic collision rate of about 60 % 
of the rate of the initial clouds. Based on this we con- 
clude that it should be possible to cool the mixed cloud 
to quantum degeneracy using evaporative cooling. 

The most critical parameter of the mixing process is 
the ratio Ii /I2 of the currents in the coils of Trap 1 and 
Trap 2. This is discussed in section fVAl There are sev- 
eral other parameters which can be varied in order to 
optimize the process. Here we look at the trap depth 
fsection [VB[) and the initial temperature (|V C[) . 



A. Dependence on trap current ratio 

Keeping I2 constant at 265 A which corresponds to a 
radial magnetic field gradient of 53 G/cm and changing 
the current ratio by varying Ii , we get the experimental 
results shown in Fig. [51 The corresponding results of the 
theoretical simulations with the same parameters lead 
to the results in Fig. [HI In these figures we plot (a) 



the number of atoms staying trapped during the process, 
(b) the final temperature, and (c) the final phase space 
density. In each of these plots the results are given for 
atoms in both traps and with atoms in only one or the 
other of the two traps. Thereby it is possible to determine 
which of the two clouds is lost or heated most during the 
process. It should be stressed that the uncertainty of 
the phase space density value is relatively large, because 
it decreases rapidly with temperature (see Eq. [7]). A 
change in phase space density by a factor of 2 corresponds 
only to a relative temperature change of 17 %. 

Both the experimental and theoretical plots show 
that for large current ratios the atoms in Trap 1 stay 
trapped during the mixing process without loss or heat- 
ing, whereas a small current ratio favors the atoms in 
Trap 2. This confirms what one could expect from the 
isosurface plots in Fig. [6l Experimentally we see that 
when we have atoms in Trap 1 only and use a current 
ratio less than 1.06, we get a diffuse cloud which is very 
elongated in the horizontal direction and not in thermal 
equilibrium. Therefore we cannot attribute a tempera- 
ture to these clouds [5l|. This is the reason why the 
"Trap 1" curve stops in the temperature and phase space 
density plots. The same situation applies when we have 
atoms only in Trap 2 at a current ratio exceeding 1.14. 
The number of atoms in these cases is determined from 
a sum of pixel values after 2 ms time-of-flight. 

Fortunately it is possible to choose a current ratio 
which constitutes a compromise where the heating of the 
mixed cloud as a whole is small and almost no atoms are 
lost. Both experimentally and theoretically there is a well 
defined current ratio which at the same time minimizes 
the temperature and maximizes the number of atoms of 
the mixed cloud. We will refer to this current ratio as the 
optimal current ratio. Experimentally it is 1.125 while in 
the theoretical calculations it is 1.215. We believe that 
the 8 % difference between these two values can be at- 
tributed to uncertainties in the calculation of the mag- 
netic field from the coils. 

The qualitative agreement between the measurements 
and the theoretical simulations is striking and also, bear- 
ing in mind that the optimal current ratio is slightly dif- 
ferent, even the quantitative agreement is not bad. Sim- 
ilar to the experimental results, the theoretical simula- 
tions give a final phase space density of about 27 % of 
the initial phase space density of the initial clouds at the 
optimal current ratio. 

But there are also differences. The simulations re- 
sult in two local maxima of the final phase space den- 
sity whereas experimentally only one clear maximum is 
observed. Also, theoretically at the optimal current ra- 
tio the phase space density of the atoms from each of 
the two traps is almost equal. Experimentally the phase 
space density is significantly larger when we have atoms 
only in Trap 1 than when we have atoms only in Trap 2 
at the optimal current ratio. Taking the simplicity of the 
theoretical model into account one cannot expect perfect 
quantitative agreement. We recall that the model does 
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Figure 8: (Color online) Experimental results of varying the 
ratio /1//2 of the currents in the coils of Trap 1 and Trap 
2, respectively {I2 is kept constant while /i is varied), (a) 
shows the number of atoms staying trapped during the process 
(normalized to the number of atoms in each of the initial 
clouds which is 1.7 x 10*), (b) shows the final temperature 
and (c) shows the final phase space density (normalized to the 
value 1.1 X lO"'^ which is the phase space density of an initial 
cloud with 1.7 X 10* atoms and a temperature of 195 /iK). All 
three plots contain three curves: One curve showing results 
with atoms in both traps and two curves showing results with 
atoms in only one or the other of the traps. The dotted lines 
between the points are drawn to guide the eye. 
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Figure 9: (Color online) Theoretical results of varying the 
current ratio. The plots show the same quantities as the plot 
in Fig. [8] Atoms which escape the trap potential or hit the 
chamber wall during the mixing process are removed before 
calculating the temperatures and phase space densities. The 
number of atoms and the phase space density are normalized 
to the (identical) values of each of the initial clouds. 



not include the dynamics before and after the mixing 
process and it ignores coUisions and interactions. 



B. Dependence on trap depth 

To the extent that gravity can be neglected, changing 
the trap depth only changes the magnitude and not the 
geometrical shape of the potential. Therefore this pa- 
rameter is not expected to be very critical as long as the 
traps are deep enough that the atoms do not escape. This 
is confirmed by our measurements. First of all there is no 
measurable change in the optimal current ratio when we 
change the trap depth. For small trap depths the final 
temperature has no well defined minimum as we vary the 
current ratio, but the number of atoms peaks at a cur- 
rent ratio which is independent of trap depth. Therefore 
it makes sense to vary the trap depth keeping a constant 
current ratio. 

The number of atoms remaining trapped after mixing 
versus the trap depth is shown in Fig. 1101 where the ex- 
perimental and theoretical results are compared. They 
agree very well. For sufficiently deep traps there is essen- 
tially no loss of atoms. When the trap gradient becomes 
smaller than a critical value which is about 60 % of the 
value used above (which corresponds to a current of 265 
A in Trap 2 giving an initial radial field gradient of 53 
G/cm and a barrier height of 154 G - see Table |T| then 
trap losses start to set in. At a trap depth of 40 % of 
the value used above only 14 % of the atoms are left. We 
have also measured how the temperature of the mixed 
cloud depends on the trap depth. In the region with es- 
sentially no loss of atoms the temperature variations are 
less than 10 %. It should be kept in mind that these re- 
sults are obtained for the |F = 1, m-p = — 1 > hyperfine 
state which has a magnetic moment of /iB/2. Using the 
\F = 2, m-p = 2 > state, which has a magnetic moment 
of /iB , the same potential energy gradient can be achieved 
with a field gradient which is only half as large. 



C. Theoretical dependence on the initial 
temperature 

To check how generally applicable the mixing proce- 
dure is, it is also interesting to look at the dependence 
on the initial temperature of the clouds which are mixed. 
We have studied this in the theoretical simulations using 
the radial field gradient 53 G/cm in Trap 2. In all the 
calculations we have used the current ratio 1.215 which 
is theoretically optimal at low temperatures (as long as 
there is essentially no loss of atoms). For higher tem- 
peratures the optimal current ratio is less well defined, 
because there is no clear minimum of the final tempera- 
ture. This is because there is a significant loss of atoms 
and the most energetic atoms have the largest probabil- 
ity of escaping. The fraction of atoms remaining trapped 
after the mixing process is plotted in Fig. [11] and the 
final temperature is plotted in Fig. [T^ 

The figures show that up to an initial temperature of 
200 /iK all atoms stay trapped whereas for higher tem- 
peratures the loss gradually increases until at 1 mK 85 
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Figure 10: (Color online) The number of atoms staying 
trapped during the mixing process relative to the initial num- 
ber of atoms in each trap for different trap gradients. The 
measurements are normalized to 1.81x10® atoms. The cur- 
rent ratio is I1/I2 ~ 1.125 for the experimental results and 
1.215 for the theoretical results. As long as the initial radial 
field gradient of Trap 2 is larger than 35 G/cm (correspond- 
ing to an initial barrier height of 92 G in Trap 2) the loss of 
atoms is only about 5 %. When the field gradient becomes 
smaller than this value a significant loss of atoms sets in. The 
experimental and theoretical results agree very well. 



% of the atoms are lost. On the other hand one does 
not gain very much from having an even colder cloud 
initially since the absolute temperature increase is about 
100 /xK and roughly independent of the initial tempera- 
ture up to an initial temperature of about 500 ^K. For 
even higher temperatures the final temperature drops be- 
low the initial temperature, because the most energetic 
atoms escape the trapping potential or hit the vacuum 
chamber wall. 



VI. OUTLOOK: DIFFERENT SPECIES 

Although of limited practical interest, mixing of two 
Rb clouds is an excellent way to study the process of 
mixing atomic clouds by merging two quadrupole traps. 
The perspective is to use the process to mix two different 
species. We hope that a cold mixture produced in this 
way can be a good starting point for making a quantum 
gas mixture. It is our intention to mix *^Rb and ^Li or 
^Li and do evaporative cooling on Rb to get both species 
into the quantum regime. 

For alkali atoms it is normally preferable to have both 
species in the stretched hyperflne state with gprnp = 1. 
Besides having the largest possible magnetic moment, the 
trapped atoms in such a mixture are least affected by 
loss due to spin-exchange collisions. Thus, assuming that 
both species are in a state with gprnp = 1, the only 
difference between the two species affecting the dynamics 
is their initial temperature (which is typically the same 
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Figure 11: (Color online) Results of simulations showing the 
fraction of atoms staying trapped during the mixing process 
for different initial temperatures. For the trap gradient used 
(53 G/cm in trap 2) there is essentially no loss of atoms until 
the initial temperature rises above about 200 /iK. At a tem- 
perature of 1000 fiK only 15 % of the atoms stay trapped. The 
current ratio is 1.215. Each point corresponds to a simulation 
using 2x2000 atoms initially. 
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Figure 12: (Color online) Theoretical dependence of the tem- 
perature after the mixing process on the initial temperature 
at a current ratio of 1.215 which is optimal for low temper- 
atures. Up to an initial temperature of about 500 /iK the 
absolute temperature increase is roughly independent of the 
initial temperature. For very low initial temperatures, the fi- 
nal temperature converges towards 150 fiK for a mixture with 
equal numbers of atoms in the two traps. At this stage the 
properties of the mixed cloud are completely determined by 
the mixing process and independent of the initial tempera- 
ture. Apparently for the trap gradients used the process is 
also feasible for higher temperatures, but then the most ener- 
getic atoms escape or hit the vacuum chamber wall. The thin 
black line shows the initial temperature for comparison. 



within a factor of two or so) and their mass which enters 
via Newton's 2"^^ law and via the initial velocity which 
has the mean value v = y^SkT/mn. 

Using our theoretical model which simulates the es- 
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Figure 13: (Color online) Simulation results of mixing of two 
different species. The plot shows the final phase space density 
of each of the clouds resulting from mixing of Rb in Trap 1 
(•) with respectively Li (o), K (t) and Cs (x) in Trap 2. The 
results for Na resembles the results for K very much. The 
temperature of the initial clouds at a radial field gradient of 
42 G/cm is 250 /iK except for Li for which it is 500 /iK. 
The phase space densities are normalized to the values of the 
initial clouds. Each point corresponds to a simulation using 
2x5000 atoms initially. 



sential features of the Rb-Rb mixing very well, we can 
make predictions on the outcome of an experiment where 
two different species are mixed and the results are very 
promising. We have simulated mixing of Rb in Trap 1 
with Li, Na, K and Cs in Trap 2. The radial field gradient 
of Trap 2 was 53 G/cm as above such that the magnetic 
trapping potential is twice as large as the potential for 
the \F = 1, mp = — 1 > ^^Rb atoms used above. We 
have chosen the initial temperature 250 /iK at a radial 
trap gradient of 42 G/cm except for Li for which 500 /iK 
is more realistic (52j . 

For current ratios in the range 1-1.45 the loss of atoms 
is below 5 % except for the Li cloud for which the loss 
is 8-12 %. Fig. [T3] shows the final phase space density 
relative to the initial phase space density of each of the 
clouds. Just as for Rb-Rb mixing the current ratio can 
be adjusted to favor one or the other of the two atomic 
clouds and there is a region of current ratios where the 
phase space densities of both species are retained at a 



reasonable level. 

Since evaporative cooling is notoriously a very lossy 
process, if one of the species is to be sympathetically 
cooled by the other, one typically wants to have many 
atoms acting as the cooling agent and fewer atoms to 
be sympathetically cooled (and it does not matter if the 
temperature of these few atoms is larger in the beginning 
as long as their contribution to the total energy is small). 
Therefore one might benefit from choosing a current ratio 
which minimizes the loss and heating of the cooling agent 
while keeping the loss and temperature increase of the 
other species at an acceptable level. 

VII. CONCLUSION 

In conclusion, we have demonstrated a new method to 
make mixtures of magnetically trapped, ultracold atoms 
based on merging of simple magnetic quadrupole traps 
mounted on industrial translation stages. This method 
should be generally applicable provided that the trap gra- 
dients and the diameter of the vacuum chamber tube in 
the mixing region are large enough. 

Our experimental results show that two clouds of Rb 
atoms can be mixed with only about 5 % loss of atoms, 
and with a final phase space density which is 24 % of the 
initial one. This loss of phase space density is so small 
that it should be possible to cool the mixed cloud down 
to quantum degeneracy. 

The experimental results are in good agreement with 
simple classical simulations of the mixing process. Ex- 
tending these simulations to mixtures of different ele- 
ments, we get promising results for mixing of different 
species and we believe that the method will become a 
useful new way of producing mixtures of ultracold atoms. 
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